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1. Introduction 



The construction and study of discrete Painleve equations has been a topic of research 
interest for almost two decades, [28], [131 ESI [21] . Reviews of the subject may be found in 
[TUdE]. The subject has culminated in the classification by H. Sakai of discrete as well 
continuous Painleve equations based on the algebraic geometry of the corresponding 
rational surfaces associated with the spaces of initial conditions [37]. As a byproduct of 
the latter treatment, a "mistress" discrete Painleve equation with elliptic dependence 
on the independent variable was discovered. 

In the history of the Painleve program, after the classification results for second- 
order first-degree equations, Painleve's students, Chazy and Gamier, [6j [7J [To] , 
investigated the classification of second-order second-degree equations and third-order 
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first-degree equations. The classification of the second- degree class was completed by 
Cosgrove in recent years, [HIS]. A partial classification for the third-order case was also 
obtained by the aforementioned authors. The work of Bureau, [Hill], is also important 
in this respect. No classification results exist for the analogous discrete case and hardly 
any examples of second-order second-degree difference equations exist to date, with the 
notable exception of an (additive difference) equation given by Estevez and Clarkson 



A key result of this letter is a second-order second-degree equation, which can be 
considered as a g-analogue of an equation in the Chazy- Cosgrove class, together with its 
Lax pair (i.e., isomonodromic g-difference problem). This new equation contains four 
free parameters, which suggests that it could be a g-difference analogue of the second- 
order second-degree differential equation that is a counterpart of the sixth Painleve 
equation. There are several forms of a second-order second-degree equation related to 
the sixth Painleve equation that have appeared in the literature, notably one derived by 
Fokas and Ablowitz [12] and another appearing in the work of Okamoto [32J. Difference 
analogues of the Fokas-Ablowitz equation have been provided by Grammaticos and 
Ramani, [HI], but these difference equations were all of first-degree. It has been argued 
by these authors that equations that are second-degree in the highest iterate cannot be 
viewed as "integrable" , however, for the new equation we establish integrability through 
a Lax pair in the form of an isomonodromic deformation system. Furthermore we show 
that the equation arises as a similarity reduction from an integrable partial g-difference 
equation. Through the same procedure we also construct higher-order second-degree 
equations, which form a hierarchy associated with the new equation. 

2. g-Difference Similarity Reduction 

Lattice equations of KdV-type were introduced and studied over the last three decades 
[T9| [29] . see [27] for a review. These lattice equations can be formulated as partial 
difference equations on a lattice with step sizes that enter as parameters of the equation. 
Conventionally we think of these parameters as fixed constants. However, in agreement 
with the integrability of these equations, there exists the freedom to take the parameters 
as functions of the local lattice coordinate in each corresponding direction. In this paper 
we consider the case when the parameters depend exponentially with base q on the lattice 
coordinates. 

We work in a space J of functions / of arbitrarily many variables Oj (i — 1, . . . , M 
for any M) on which we define the g-shift operations 



For u,v,z G 9 r , we consider the following systems of nonlinear partial g-difference 
equations: 



urn- 



£Tif(ai, ...,a N ):= f(ai 



, . . . , gaj, . . . , ajvf) • 





(2.2) 
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and 

Oj (z qT%Z) ( qTjZ qTi qTjZ*) 

= aj(z- qTjz) ( qTiZ - ^ qTjz) , (2.3) 

where i,j = 1,...,M. Each of these systems, (12.11) to (12.31) . represents a multi- 
dimensionally consistent family of partial difference equations, in the sense of [3U [T], 
which implies that they constitute holonomic systems of nonlinear partial g-difference 
equations. Another way to formulate this property is through an underlying linear 
system which takes the form 

q Tr 1 <j> = M i (k)(f>, (2.4) 

where <f> = <fi(k; {%}) is a two-component vector-valued function and by consistency, 
q T~ q Tf l (f) = q Tf l q Tf l (p, leads to the set of Lax equations (for each pair of indices 

( q T^M 3 )M t = ( q T^ x Mi)Mj . (2.5) 

We will consider three different cases, associated respectively with equations ( 12. ip - 
(12. 3p . To avoid proliferation of symbols we use the same symbol Mi(k) for each of the 
respective Lax matrices. For specific choices of the matrices Mi the Lax equations (12.51) 
lead to the nonlinear equations given above. In the case of the g-lattice KdV, ( 12. ip . the 
Lax matrices Mi are given by 

In the case of the g-lattice mKdV, (12.21) . the Lax matrices Mi are given by 

M,(t, {«,}) = -L, ( a '%" )/V • *' v ) . (2.7) 

Finally, in the case of the g-lattice SKdV, (12. 3p . the Lax matrices Mi are given by 

M((M<lj}) « ( * ' W4H'-fT*')- 1 ). (2.8) 

Wj — K \ Z — q L i Z ^ 1 / 

These Lax matrices are straightforward generalizations of those with constant lattice 
parameters given in e.g. [30] . 

We mention that the solutions of the equations (12. ip to (I2.3P are related through 
discrete Miura type relations, namely 

ai (z- q T- l z)=v( q Tr 1 v) , (2.9a) 
s = (a* - qT^u) v - ^ qT~ l v , (2.96) 
qT^s = a t v - (ai + u) gT^v , (2.9c) 

where s G J is an auxiliary dependent variable. From these relations, the partial q- 
difference equations ( 12.11) to ( 12.31) can be derived by eliminating s. 

Similarity reductions of lattice equations have been considered in [281 E5J [31j [26J, [30] 
where it was shown that scaling invariance of the solution can be implemented through 
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additional compatible constraints on the lattice equations. In the present case of (12.11) 
to (12.31) these constraints adopt the following form [11] 

u({q- N ai}) = q~ N ^ ^ > -^— ] (2.10a) 

KVH lJJ H l + A(g w - l)(-l)Si 9l °s a * Vl JJJ v ; 

1 _ \( n N _ 1U_1 ^Ei 'log »i 

«({?-^ai}) = i + u « a i})> ( 2 ' 106 ) 

where A and /i are constant parameters of the reduction and N £ N represents a 
"periodicity freedom" . The notation 9 log x refers to the logarithm of x with base q. 

In order to compute the corresponding isomonodromic deformation problems 
associated with the similarity reductions we have the following constraints on the vector 
function of the Lax pairs. In the case of (12.11) we have 

<t>(q N k;{a 3 })= (2.11) 
(l + A(^-l)(-lp"^) ,0 u 
-2A^ (£, a,) (-lp , ^(l-A(^-l)(-l)^-) a ^>- 

In the case of (J22D 

(t>(q N k;{ aj })= (2.12) 

(l + A(g JV -l)(-l)^ 91 °s^) 
q- N (l+n(q N -1)) 

In the case of (12.31) 

4>(q"k;{a j })= (2.13) 

(l-M^-i)) o 

q- N (l+fi(q N -1)) 

The similarity constraints, (12. lip to (12.131) . in conjunction with the discrete linear 
equations (12. 6p to (12.81) can be used to derive corresponding g-isomonodromic 
deformation problems. That is, (12. lip to (I2.13P lead to g-difference equations in the 
spectral variable k, hence together with the lattice equation Lax pairs we obtain q- 
isomonodromic deformation problems for the corresponding reductions. 



#;{r v %}). 



0(*;{g-%}) 



Remarks: 

(i) The similarity constraints above were obtained through an approach based on 
Jackson-type integrals, the details of which will be presented elsewhere [TT|. By 
construction, these constraints are compatible with the underlying lattice equations, 
which can be checked a posteriori by an explicit calculation, presented in the 
appendix. 

(ii) In this approach, the dynamics in terms of the variables ai appear through 
appropriately chosen g-analogues of exponential functions, whereas the relevant 
Jackson integrals exhibit an invariance through scaling by factors q N . 
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(iii) The parameters A and fi arise in this setting through boundary contributions in a 
manner analogous to the derivation in [30] . 

In the remainder of this letter our aim is to implement the similarity constraint to 
obtain explicit reductions to ordinary g-difference equations. For simplicity we consider 
only the reduction of the g-mKdV equation ( 12. 21) . leaving considerations of the q-KdV 
and g-SKdV to a future publication [II] . There are two possible scenarios to derive 
similarity reductions of the lattice equations using the constraint (12.10&j) . 

"Periodic" similarity reduction: By fixing M = 2 and allowing N to vary, we select 
two lattice directions, say the variables a\ and 02, and consider similarity reductions 
with different values of N. This is a g-variant of the periodic staircase type reduction 
of partial difference equations on the two-dimensional lattice. For instance, with N = 2 
the reduction is a second-order first-degree g-Painleve equation. Increasing iV leads to 
g-difference Painleve type equations of increasing order. However, we will not pursue 
this route here but leave it to a subsequent publication [II]. These reductions are 
reminiscent of the work [181 ESI 02] • 

Multi-variable similarity reduction: The similarity constraints provide the mechanism 
to couple together two or more lattice directions. By considering the case N = 1 we 
implement the similarity constraints on an extended lattice of three or more dimensions 
in order to obtain coupled ordinary g-difference equations, in a way that is reminiscent 
of the approach of [31]. This is considered in the next section. 

We have not considered the more general case of arbitrary M, N e N, which we will 
postpone to a future publication [IT] . 

3. Multi-variable similarity reduction 

In this section we consider explicitly the M — 1, 2 and 3 cases for N = 1. 

For simplicity we shall in what follows denote the coefficient in ( 12.1051) as 7, i.e. 



In contrast to the usual difference case which was explored in [31] where in the case 
of two variables we obtain a nontrivial OAE as a reduction, in the g-difference case we 
have to consider at least three independent variables to obtain a nontrivial system of 
OAEs reduction. 

In [31] the compatibility between the similarity constraint and the lattice system 
was established and led to a system of higher order difference equations in the reduction, 
namely equations which were on the level of the first Gamier system. In contrast to 
the q = 1 work, the 3D similarity constraint here is somewhat simpler and leads to a 
second-order equation (which is of second-degree, and is a principal result of this letter). 




v({q ai}) = 7^(R}) 



(3.1) 



where 7 alternates between two values, i.e., ^Tf r j = 7. 
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Two-variable case 

Let us now select among the collection of variables {aj} two specific ones which for 
simplicity we will call a and b. Denote the g-shifts in these variables by an over-tilde 
and an over-hat respectively. Equation (12.21) may now be written 

bvv + avv = avv + bvv, (3.2) 

where the over-tilde ~ refers to the g-translation a ^ qa and the over-hat ^ refers to the 
g-translation b^> qb (so if v = v(a, b),v = v(qa, b), v = v(q~~ l a, b), v = v(a, bq), . . . ). 

Equation (13. 1 p gives the constraint v = yv to impose on (13.21) (where 7 = 7 = 7). 
This leads to the linear first-order (in that it is a two point) ordinary difference equation 

v = Cv, (3.3) 

where C = 7(07 + b)/(a + 67). In the appendix the consistency between the lattice 
equation (13. 2p and the constraint (13.31) is shown by direct computation. 

Three-variable case 

Take three copies of the lattice mKdV equation with ai = a, a 2 = b, a 3 = c, 

bvv + avv = avv + bvv, (3.4a) 
cvv + avv = avv + cuv, (3.46) 
cvv + bvv = bvv + cvv, (3.4 c) 

where c = qc, together with the constraint 

v(g _1 a, q~ l b, q~ l c) — jv(a,b,c). (3.5) 

(The similarity constraint is shown by a direct computation to be compatible with the 
multidimensionally consistent system of mKdV lattice equations in the appendix.) We 
now proceed to derive the reduced system which leads to a (higher-degree) ordinary 
g-difference equation in terms of one selected independent variable, say the variable 
a. The remaining variables b and c will play the role of parameters in the reduced 
equation. Thus, we can derive the following system of two coupled OAEs for v(a, b, c) 
and w(a, b, c) = v(a, b, q~ x c): 

arjv — bw 

jv = w — - , (3.6a) 

a w —b'jv 

aw — cv 

w = v , (3.66) 

a v —cw 



where 7 = 7. We consider the system (13.6 ap and (13.6fe|) to constitute a g-Painleve system 
with four free parameters. 
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The system (13.6 ap and (13.661) can be reduced to a second-order second-degree 
ordinary difference equation as follows. Introduce the variables 



X=- , V 
w 

then from (13.6 ap we obtain 
7- = jVXW 



v 

V 



W 



w 
w 



w 



whereas from (13.66P we get 



a>yX + 6 
kyX + a 



q~ l cXV 



(3.7) 



(3.8) 



(3.9) 



X q^c/X + aV 

using also the definitions (13. 7p . Thus, we obtain a quadratic equation for V in terms of 
X and X and hence also we have W in terms of X and X. Inserting these into (13.81) we 
obtain a second-order algebraic equation for X. Alternatively, avoiding the emergence 
of square roots, the following second-order second-degree equation for X may be derived 

N on 2 

2 ~ v , a~{X + b 



rxx 



bjX 



,c 2 1 



/ a'yX + b 



1 a 2 X \lrfX 



+ a 

g -1 W(i 



XX) H 



XX) + q-\l-XX) 



a>yX + b 
byX + a 



XX) 



a~{X + b 



(3.10) 



+ a 

We consider this second-degree equation to be one of the main results of this letter. 

We now proceed to present the Lax pair for the g-Painleve system (13.6 a-j) and 
( 13.66P and the second-order second-degree equation ( 13 . 1 Of) . The Lax pair is formed by 
considering the compatibility of two paths on the lattice: along a 'period' then in the a 
direction and evolving in the a direction then along a 'period'. Using (12. 12ft the evolution 
along a period is converted into a dilation of the spectral parameter, k, by q. The result 
is the following isomonodromic g-difference system for the vector <j)(k; a) which using 
the results of section [2] yields 



where 



(j)(k; q 1 a) 
4>{qk; a) 



M(k;a) 



M(k;a)(f)(k;a) 
L(k; a)(f)(k; a) 



av/v k 2 /v 



(3.11a) 
(3.116) 



(3.12a) 



and 



L(k; a) 



ajv/v k 2 /v 
q~ l ^v q~~ l a 



byv/w k 2 /w 
b 



cw/v k 2 /v 



w 



(3.126) 
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where we have suppressed the dependence on the variables b and c (which now play the 
role of parameters) and omitted the unnecessary prefactors (6 — k)~ l and (c — k)~ l , as 
well as an over factor g _1 (l + /i(g — 1)). 

The consistency condition obtained from the two ways of expressing <p(qk; g -1 a) in 
terms of <j)(k; a) is formed by the Lax equation 

L(k; q~ 1 a)M{k; a) = M(qk; a)L(k; a) . (3.13) 

A gauge transformation can be obtained expressing the Lax matrices in terms of the 
variables introduced in (13.71) . Setting 

M(*;«) = _^ 1 (3.14a) 

, , 1 / ^(ab^X + k 2 ) k 2 ( ai X + b)/V \ ( c/X k 7 
[ ' flj ~ a-k\ q-^Via + bqX) q-\ab + k^X) J \ 1/X 

the Lax equations (I3.13P (replacing L and M by £ and M respectively) yield a set of 
relations equivalent to the following two equations: 

~m = 22*±* , (3.i5) 

a + tryX 

aV 2 + q- 1 c(^-X\V-a^ = , (3.16) 
using also 7 = 7. This set follows directly from (13. 8p and (13.91) . Thus A3. 14oj) and (13.1461) 
form a g-isomonodromic Lax pair for the second-degree equation (13.101) . 

Four-variable case: 

Suppose we have 4 variables a iy i = 1,...,4. Select a — a\ to be the independent 
variable after reduction. Introduce the dependent variables Wj_ 2 = qTf v, j = 3, 4. 
Then directly from the g-lattice mKdV equation (12. 2p we have the set of equations 
awj — Oj +2 v 

Wj = v , j — 1,2 , (3-17) 

a v —aj + 2Wj 

where as before the tilde denotes a g-shift in the variable a. At the same time the 
multiply shifted object ^Tg -1 ^T^ 1 v can be expressed in a unique way (due to the CAC 

property) in terms of v and ^f 1 ^ = iWj_2, j = 3,4, by iterating the relevant copies of 

the g-lattice mKdV equation in the variables aj, j 7^ 2, leading to an expression of the 
form ^r 3 _1 gT^ 1 v =: F(v, Wi, W2), where F is easily obtained explicitly. Imposing the 

similarity constraint (13.11) we obtain 7 ^T 2 v = F(v, W\, w 2 ) and inserting this expression 

into the g-lattice mKdV ( 12. 2p with % — l,j = 2 we obtain 

a 2 a 3 w 2 - a A wi \ ( J\ 2 2n ~ / Q1fi x 

a H a 2 7 r [v, W\,w 2 ) — av \ = (a 2 — a )vv . (3.18) 

7 a 3 W! - a 4 w 2 J V / 



, (3.146) 
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(a 2 - a 2 )77^(a 3 u;i - a 4 w 2 ) 



a(a\ — a 2 ) v +a 3 (a 4 — a 2 )wi + a 4 (a 2 — a\)w 2 



[(a 2 a 3 - ^aa A )w 2 + (7aa 3 - a 2 a 4 )u>i] 



a(a 2 — a 2 A )(a 2 w\w 2 — 'javv) 



+ ( a 2 a 4 (a — a 3 ) v—aa 3 (a 4 — a )jv J W\ + ( a 2 a 3 (al 



v j w 2 , 
(3.19a) 



and this is supplemented by the two equations 

avwi + a-sWiWi = a 3 vv +av w\ , 

avw 2 + a i w 2 w 2 = a 4 v v +a v w 2 , 



(3.196) 
(3.19c) 



which is equivalent to a five-point (fourth-order) g-difference equation in terms of v 
alone, containing five free parameters: a 2 , a 3 , a 4 , A and /i (inside 7 and 7). This would 
be an algebraic equation, so we proceed as follows in order to derive a higher-degree 
g-difference system. Introduce the variables 

X i = - , Wi = ^ , i = l,2, 

Wi Wi 

while retaining the variable V = v/v as before. By definition we have 

V X 



(3.20) 



z = 1,2, 

Wi Xi ' 

and from H3. 19 6|) . (13. 19 eft we obtain 

qa + a i+2 VXi VX { 



Wi 

qaV + a i+2 /Xi 
whilst from ( 13. 19 aft we get 



Xi 



% = 1,2 



(3.21) 



(3.22) 



(a 2 2 - a 2 ) 7 7^ 



a 3 



X 1 X, 




a 4 — a a — c 

+ a3 x l + a4 ^X 



o-2 a 3 — jaa^ jaa 3 — a 2 a^ 



Xo 



X 1 



f 2 2\/ fl 2 



-lay, 



,2 „2 



+ a 2 a 4 



Ob ^3 / 2 2 V t / 1 ^" I ^4. ^ - ■ /O o\tt-\ 1 



- aa 3 (a 4 - a )7V J — + | a 2 a 3 - 



— 7aa 4 (a — a 3 )V 



From (I3.22p we obtain the set of quadratic equations for V 

i—V + a i+2 -=r 
X: VX 



gai^K 2 + a i+2 (i -X I F-ga = , i=l,2 



(3.24) 



X 2 
(3.23) 
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from which by eliminating V we obtain 



a 3 (l - XxX^Xs - 04(1 - X 2 X 2 )X 1 
= q 2 a 2 (X 1 X 2 -X 2 X 1 ) 2 . 
Furthermore, solving V from the quadratic system as 

X 2 Xi — X\X 2 



a 3 (l - X 1 X 1 )X 2 - a 4 (l - X 2 X 2 )X 1 



V 



qa- 



(3.25) 



(3.26) 



a 3 (l - XxX^Xa - a 4 (l - X 2 X 2 )X 1 
and inserting this into (13. 23ft we obtain a second-order equation in both X%, X 2 coupled 
to the equation (13. 25j) which is first order in both X±, X 2 . It is this coupled system of two 
equations in X±, X 2 which forms our higher order generalisation of ( 13.101) . The system 
of (13.231) and (13.251) with (13.261) constitutes a third-order system with five parameters. 

The derivation of the Lax pair follows the same approach as that for the three- 
variable case (with an extra factor in £ due to the additional lattice direction). We 
omit details here, which we intend to publish in the future pTj . 



Beyond the four-variable case: 

It is straight-forward to give the form of the full hierarchy, however due to lack of space 
we postpone this until a later publication [11] . 



4. Conclusion and discussion 



In this letter we have presented the results of a scheme to derive partial g-difference 
equations of KdV type and consistent symmetries of the equations and demonstrated 
how it can be implemented. Lax matrices follow from this approach. A notable result is 
the derivation of the higher-degree equation (13.101) . showing that the scheme presented 
here allows for the derivation of new results within the field of discrete integrable 
systems. 

The first-, second- and third-order members of the N — 1 hierarchy have been 
presented. The scheme continues to give successively higher-order equations by 
considering successively higher dimensions of the original lattice equation. One may ask 
the natural question as to whether this gives an 'interpolating' hierarchy which, contrary 
to the usual cases, increases the order and number of parameters of the equations by 
one in each step, rather than a two step increase. A further natural question connected 
with this hierarchy is its relation to the g-Garnier systems of Sakai [38J. 

We will present full details of the scheme from which the lattice equations (12. ip to 
( 12. 3p and their associated constraints follow in a future publication [XT]. There we will 
consider the most general case of symmetry reductions (arbitrary N £ N) of all three 
lattice equations. 

We also intend to return in a future publication to the question of limits and 
degeneracies of the equations presented in this paper. These include the q — > 1 
continuum limit, the q — > 1 discrete limit and the q — ► crystal or ultradiscrete limit. 
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5. Appendix 

In this appendix we present the result of explicit calculations showing the consistency 
of the lattice equations and constraints. 



Two-variable consistency 

We shall check the consistency between the lattice equation f l3.2j) and the constraint 
v = yv by direct computation. This computation is illustrated in the following diagram: 



V2 
X- 



■0 "12 



V. 



Fig 1. Consistency on the 2D lattice. 

Assuming the values vo, v\ as indicated in Fig 1 are given, we compute successively v 12,1*2 
etc., where the subscripts refer to the shifts in the lattice variables a,b respectively, as 
is evident from Fig 1. Points other than vq and V\ are computed using either the lattice 
equation (indicated by x) or by using the similarity constraint (indicated by O)- The 
value v_i_2 is the first point which can be calculated in two different ways (hence 
indicated in the diagram by <8>). Without making any particular assumptions on how 7 
depends on a and b, a straightforward calculation shows that the two ways of computing 
V-1-2 are indeed the same, for any choice of initial data vo and v%, provided that 7 obeys 
the relation 



a + 67 
6 + 07 



bj 



07 



7 
7 



(5.1) 
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A particular solution of this relation is 

7 = 7 7 = 7) ( 5 - 2 ) 

and hence 7 = 7 implying that 7 is an alternating "constant" which is in accordance 
with the value given in (12.10ftj) . The reduced equation in this case is (13. 3p . which can 
be readily integrated. 

More generally, equation (15.11) can be resolved by setting 

a + h v v . , 

7- = - , 7 = - , (5.3) 

+ 07 v v 

leading to the consequence that v has to solve the q- lattice mKdV (13.21) . In principle we 
could take for v any solution of the reduced equation (13. 3p and use this to parametrise 
the reduced equation for v via the relations ( 15. 31) . In any event, we see that the two- 
variable case does not lead to interesting nonlinear equations. 



Three-variable consistency 

In this case the consistency diagram is as follows: 



^0 Vi 




V-l-3 V-3 

Fig 2. Consistency on the 3D lattice. 

A similar notation as the previous case is used as is evident from Fig 2. The initial data 
vq, vi and v^ are given, and the indicated values on the vertices are computed either by 
using one of the lattice equations ( 13. 4 aft to ( 13. 4 eft or the similarity constraint ( 13.51) over 
the diagonal. Thus, V\p is obtained from (13.4aj) yielding 

av 2 — bv\ 

Vl,2 = V — , 

avi — ov 2 

whilst from the similarity constraint we obtain 

V-1,-3 = I2V2 , V-3 = 71,2^1,2 , V-2,-3 = llVl , 

assuming that 7 shifts along the lattice, indicated by the indices, and finally the value 
of t>_i 5 _2,_3 can be computed in two different ways, leading to 

av-2,-3 - &y-i,-3 ajtVi - 672^2 av 2 - bvi 

u-1,-2,-3 = 7o^o = r v -3 = r 7i,2^o 7— , 

av-i 3 - ov-2-3 072^2 - (ryivi avi - bv 2 
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leading quadratic identity in Vi and v 2 . Assuming that the latter must hold identically, 
and thus setting all coefficients equal to zero, we obtain the following conditions on 7: 

71,2,3 = 7i = 72 = 73 , 

from which we conclude that 7 is an alternating "constant" , for instance 

7 = a p ■ ' (a, p constants) (5.4) 

and this leads to the conditions from which it is easily deduced that the form (13. ip of 7 
satisfies these conditions. 
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